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ABSTRACT 


A  FORTRAN  program  has  been  developed  to  allow  for  the 
use  of  probabilistic  design  methods  in  the  numerical  optimi- 
zation process.   The  program  was  written  as  a  set  of  sub- 
routines for  COPES  (Control  Program  For  Engineering  Synthesis) . 
COPES  maximizes  or  minimizes  a  numerically  defined  objective 
function  subject  to  a  set  of  inequality  constraints  using 
the  optimization  program  CONMIN  (A  Fortran  Program  for  Con- 
strained Function  Minimization) .   The  program  developed  here 
allows  for  the  use  of  both  the  normal  and  lognormal  distribu- 
tion models.   Design  examples  are  presented  to  demonstrate 
the  program  capabilities.   User  instructions  are  provided  for 
inclusion  in  the  COPES  user's  manual. 
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I.   INTRODUCTION 

In  the  traditional  deterministic  approach  in  engineering 
design  the  variables  are  treated  as  single  valued  numbers. 
These  variables ,  whether  dimensions,  material  properties, 
loads,  etc.,  are  actually  statistical  in  nature.   By  the  use 
of  safety  factors  designers  protect  against  these  variations 
by  usually  over  designing.   A  more  logical  approach  is  to 
take  into  account  the  statistical  data  known  about  each 
variable  and  design  for  a  certain  reliability.   By  consider- 
ing the  statistical  nature  of  each  variable  we  should  be 
able  to  better  predict  reliability  and  performance. 

Numerical  optimization  has  proven  to  be  a  very  powerful 
tool  in  engineering  design.   Virtually  all  design  problems 
require  minimization  or  maximization  of  some  objective.   For 
the  design  to  be  acceptable,  it  must  also  satisfy  a  certain 
set  of  specified  requirements  called  constraints.   If  these 
constraints  are  specified  as  probabilities  of  failure  the 
result  would  be  a  probabilistic  design. 

The  purpose  of  this  research  was  to  test  the  applica- 
bility of  combining  probabilistic  design  concepts  with  those 
of  numerical  optimization  by:   (1)  developing  a  pilot  computer 
code  to  calculate  probabilities  of  failure,  (2)  incorporate 
this  computer  code  with  COPES/CONMIN  in  order  to  perform 
probabilistic  optimization,  and  (3)  test  using  numerous 
examples . 


11 


It  is  assumed  throughout  this  discussion  that  the  reader 
is  familiar  with  the  use  of  the  COPES/CONMIN  optimization 
program  [Ref.  1],  and  [Ref.  2]. 
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II.   NUMERICAL  OPTIMIZATION 

The  numerical  optimization  problem  considered  here  is 
stated  as  follows:   Find  the  set  of  n  design  variables 
contained  in  the  vector  X  which  will 

Minimize   F(X)  (1) 

Subject  to: 

g.  (X)  <_     0     j   =   l,m  (2) 


x£   <   x.   <   xV     i   =   l,n       (3) 
l-i-i 


The  components,  x.,  of  X  are  referred  to  as  design  variables 
which  are  changed  to  improve  the  design.   The  function  F(X) 
is  called  the  objective.   Inequality  constraints,  g.(X),  are 
the  response  limits  imposed  on  the  design.   There  are  n  de- 
sign variables  and  m  inequality  constraints.   The  lower  and 

I  u 

upper  bounds,  x.  and  x.  are  limits  imposed  on  the  design 

variables  to  insure  a  practical  result. 

For  a  deterministic  design  COPES  [Ref.  1],  determines 

the  constraints  and  objective  as  follows.   Consider  the 

design  of  a  single  bar  undergoing  uniaxial  tension.   If  we 

wish  to  minimize  the  weight,  the  objective, 


F(X)   =   pAL  (4) 
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where 

p  =  the  specific  weight; 
A  =  the  area  of  the  bar; 
L   =   the  length  of  the  bar. 

The  design  variable  is  the  area,  A.   If  the  stress  in  the 
bar  is  S  and  the  stress  limit  imposed  is  Sy  we  desire  that 
S  be  less  than  or  equal  to  Sy  so  the  constraint  becomes 

S/Sy  -  1   <   0  (5) 

In  this  investigation  a  computer  subroutine  was  developed 
to  provide  constraints  based  on  probabilities  of  failure. 

If  the  allowed  probability  of  failure  is  PF  and  the 
actual  probability  of  failure  was  calculated  to  be  Pf,  then 
the  constraint  would  be 

Pf/PF  -  1   <   0  (6) 
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III.   PROBABILISTIC  DESIGN 

A.   RANDOM  VARIABLES 

If  one  has  a  large  population  of  n  elements  for  which 
some  parameter  x.,i5=l,n  is  determined,  the  mean  value  of 
the  population  is 


ux  '    5  X  (xi>  (7> 


The  variance  is 


l   n  o 

V    =  -      I       (x.  -  y  )  (8) 

x     n  .  *•_    i   Mx 


The  standard  deviation  is 


a    =   /V  (9) 

xx 


The  standard  deviation  is  a  measure  of  the  amount  of  varia- 
bility of  the  population  data.   An  additional  measure  of 
the  variability  is  the  dimensionless  coefficient  of  variation, 


a 
C   =  —  (10) 

x     ^x 


In  engineering  design  the  coefficient  of  variation  is  fre- 
quently used. 
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B.   FUNCTIONS  OF  SEVERAL  VARIABLES 

If  Y  is  a  function  of  n  variables,  the  standard  deviation 
of  the  function  Y  can  be  approximated  as  follows: 


If   Y   =   f(X;L;  i  =  l,n) 
Assuming  that  all  x.  are  independent  random  variables 


then  according  to  [Ref.  3,  p.  59] 


1=1   1     1 


C.   STATISTICAL  MODELS 

1 .   Normal  Distribution 

The  normal  distribution  model  is  one  of  the  most 
widely  used.   The  normal  density  function  is  given  by 


where : 


fx(x)   =     ±        exp[ 2~^ — ]  (12) 

v/2tt<j         2  a 

X  X 


f  (x)   =   the  probability  density  function  of 
the  random  variable  X; 


x 


ux 


=   standard  deviation  of  the  random 
variable  X; 

=   mean  of  the  random  variable  X. 
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The  magnitude  of  the  standard  deviation  determines  the  dis- 
persion of  the  distribution.   A  large  standard  deviation  will 
result  in  a  wide  bell  curve  and  a  small  standard  deviation 
will  result  in  a  narrow  one.   The  curve  is  symmetric  about 
the  mean.   Approximately  68.26%  of  all  samples  will  lie 
within  plus  or  minus  one  standard  deviation  of  the  mean  value 
The  standardized  form  of  the  normal  density  function  [Ref.  4: 
p.  194],  can  be  expressed  as 


MY>   ■  —  exp[3L-]  (13) 

Y        /27      z 


where 


y  = 


x  '  yx 

ax 


(14) 


)    (y)   =   the  standardized  form  of  the  single 
variate  normal  density  function. 


The  single  variate  cumulative  normal  distribution  function 
is  given  by 


X 
Fx(x)   =     /   fx(x)dx  (15) 


where 


Fv(x)   =   the  probability  that  the  random  variable 
X  is  equal  to  or  less  than  the  specific 
value  x. 


X 
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Or,    using  the  expression, 


y  -  ^ 


(16) 


Equation  (15)  can  be  rewritten  as,  [Ref.  4:  p.  195] 


or 


(x-yx)/ax       2 

FX(x)   =   -=   /  exp[^-]dy 

/2tt  -oo 


x  -  y 
Fx(x)   =   $[—-£]  (17) 

x 


2 .   Lognormal  Distribution 

If  ln(X)  has  a  normal  distribution  then  X  is  said  to 
have  a  lognormal  distribution.   The  lognormal  density  func- 
tion, [Ref.  4:  p.  196],  is  given  by 


~  2 
c    i     \  11      r-(ln(x)-X)  ,  ^    n  ,,. 

fx(x)   =      -   x  exPt ^2 ];   x        ( 

/2ttg  2a 

x  x 


where 


X   =   mean  of  ln(X)  ; 


a    =   standard  deviation  of  ln(X) 

x 
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The    cumulative    lognormal   distribution    function,     [Ref.    4: 
p.    196],    is    given  by 

ln(x)-X 

1  °x  2 

FX(X)       =      T=-        I  exp[-2-]dy;    x    >    0  (19) 

/2tt    -<»  z 


or 


Fx(x)       =      *y[iHMzij  ,20) 

°x 


The   mean,    u    ,    and   the  variance   V    ,    of    X   are 
x  x 

-2 

Qx 
Ux      =      exp[X  +   -j-)  (21) 


V        =      exp[2X+ah  [expiah  -  1]  (22) 

«  xx 


Solving  for  X  and  a 
3  x 


X      =      ln(y*)    -   |  In  (Vx+u^)  (23) 


Sx      =      J-*^    +   ln[Vx+^x]  (24) 


D.       PROBABILITY    OF    FAILURE 

Consider   the    constant   area  bar    in   Figure    3.1    under   load 
P.       If   the   stress    is    the  bar    is    S    and   the    stress    limit    in 
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A 


V 


S  = 


_P 
A 


Figure  3.1.   BAR  UNDER  UNIAXIAL  LOAD 

tension  for  the  material  is  Sy,  then  the  probability  of 
failure  of  the  bar  is  the  probability  that  the  stress  S  will 
be  greater  than  the  limit  Sy.   The  method  for  determining 
the  probability  of  failure  in  the  bar  for  both  the  normal 
and  the  lognormal  formats  is  presented  below. 
1 .   Normal  Distribution 

If  both  S  and  Sy  are  random  variables  which  follow  a 
normal  distribution  then  the  probability  of  failure  is  the 
probability  the  S  is  greater  than  Sy.   The  failure  function, 
[Ref.  5:  p.  8-24] 


Z   =   (Sy  -  S) 


(25) 
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is  also  normally  distributed.  The  probability  of  failure 
is  the  probability  that  the  failure  function  is  less  than 
or  equal  to  0.0  or 

Pf   =   P(Z<_0)  (26) 

In  the  standard  form 


Pf   =   <f>  [— £]  (27) 

3  °z 


Where,  by  the  algebra  of  functions,  [Ref.  3] 


uz   =   ^Sy  "  yS  (28: 


and 


°Z   -   \/°Sy  +  °S  (29 


The  safety  index,  [Ref.  5:  p.  8-24],  is  defined  as 


wz      1 
3   =   —  =   j±-  (30 


The  probability  of  failure  is  then 


Pf   =   *g[-3]  (3D 
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or  the  area  under  the  unit  normal  density  function  from  6 
to  a.  . 


oo        .  2 

Pf      =        /      -~:  exp[^f-]d|3  (32) 

S         /2tt  ^ 


2 .   Lognormal  Distribution 

If  the  same  two  variables  S  and  Sy  follow  a  logri- 
thetic  normal  distribution  then  the  probability  of  failure 
is  expressed  by,  [Ref.  5:  p.  8-27] 

Pf   =   P(Sy^S)   =   P(ln(Sy/S)  <  0)  (33) 

Letting 

Z   =   ln(Sy/S)  (34) 


then 


Pf   =   P(Z  <  0)  (35) 


or 


Pf   =   $0[-6l  (36) 


The  safety  index,  [Ref.  5:  p.  8-29]  is 
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3   = 


ln(Sy/S) 

'ln[(l  +  C^y)  U  +  c|)  ] 


37) 


where  the  Sy  and  S  are  median  values  and 


Sy   = 


^y 


ii   +  c 


Sy 


(38) 


where 


S   = 


1  +  c 


(39. 


'Sy 


'sy 


(40) 


(41) 


3  will  be  normally  distributed  and  the  probability  of 
failure  is 


Pf   =   *, 


-ln(Sy/S) 


'ln[(l+c|  )  d+Cg)  3 


] 


(42) 


The  probability  that  S  is  greater  than  Sy  is  then  determined 
in  the  same  manner  as  the  normally  distributed  case,  by 
integrating  the  unit  normal  density  function  from  3  to  °° . 
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IV.   PROBABILISTIC  OPTIMIZATION 

The  objective  of  probabilistic  optimization  is  to  mini- 
mize or  maximize  the  mean  value  of  a  given  function  subject 
to  a  set  of  constraints  based  on  allowed  probabilities  of 
failure. 

A.   CONSTRAINTS 

The  computer  subroutine  developed  in  this  investigation, 
C0PE19,  calculates  the  constraints  as  described  in  Chapter 
II.   Given  a  constrained  variable  in  terms  of  a  desired 
probability  of  failure,  one  needs  only  to  calculate  the 
actual  probability  of  failure  to  determine  the  constraint 
value. 

Consider  the  cantilevered  beam  in  Figure  4.1.   A  proba- 
bility of  failure,  PF,  is  desired  for  stress  in  bending. 
The  mean  values  and  coefficients  of  variation  for  the  yield 
strength,  Sy,  the  length  L,  the  load  P,  and  the  dimensions  B 
and  H  are  given.   The  stress  in  bending 


SB   =   6PL  (43) 

BHZ 


The  standard  deviation  of  the  bending  stress 


/T^SB ,  2,  ,aSB,  x  2,  ,3SB^  ,  2,  ,  9SB^  .2  ,  ... 

°SB   -   W(-3P"V  +(-3lTJL)  +(V3H)  +(WV  (44) 
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Figure    4.1.       CANTILEVERED   BEAM 
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The  safety  index 


Sy-  SB 


aSy  +  aSB 


(45) 


The  calculated  probability  of  failure 


Pf  =   *g[-6]  (46) 


The  constraint  will  be 


Pf 
G   =   |f  -  1  (47) 


In  the  computer  program  developed,  C0PE19,  the  required 
partial  derivatives  of  the  constrained  variables  with 
respect  to  the  variables  in  the  expression  are  calculated 
by  finite  difference  methods. 

B.   OPTIMIZATION 

Consider  the  three  bar  truss  in  Figure  4.2.   The  geometry 
is  specified  and  there  are  two  independent  load  conditions. 
The  design  task  is  to  determine  the  areas  required  for  each 
bar  that  will  yield  a  minimum  mean  value  of  the  structure 
weight  while  keeping  the  stress  in  any  member  less  than  the 
yield  stress. 
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THREE  BAR  TRUSS 


P2 


/\ 
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Figure  4.2.   THREE  BAR  TRUSS 
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Minimize:   W   =   pH(/2Al  +  A2  +  /2A3) 
Subject  to:   Sy  compression    SIG.  .  <_  Sy  tension 
where  SIG^  is  the  stress  in  member  i  under  load  condition  j 
Given:   Material   p   =  .1  lb./cu.  in. 
Geometry   H   =   10.0  in. 
Loads  PI   =   P2   =   20,000  lb. 
Stress  limits:   Sy  compression   =   -15,000  psi 

Sy  tension   =   20,000  psi 
Maintain  symmetry  Al   =  A3 

Beginning  with  the  design  Al  =  A2  =  A3  =  1.0  sq.  in.  and 
optimizing  using  COPES/CONMIN  (deterministic)  the  following 
results  are  obtained: 

W   =   2.6  32  lb. 

Al   =   A3   ■   .7796  sq.  in. 

A2   =   .4275  sq.  in. 

The  critical  constraints  are  SIGH  and  SIG32 ,  which  are 
equal.  The  design  space  is  shown  in  Figure  4.3.  As  no 
safety  factor  was  used,  the  above  deterministic  design 
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Figure  4.3.   THREE  BAR  TRUSS  DESIGN  SPACE  (DETERMINISTIC: 


29 


corresponds  to  a  probability  of  failure  of  0.5  or  a  50% 
chance  that  SIGH  and  SIG32  will  be  greater  than  the  stress 
limit  Sy.   Utilizing  the  subroutine  COPE19  to  provide  proba- 
bilistic based  constraints  and  assuming  that  the  coefficients 
of  variation  of  all  the  variables  are  0.1,  the  following 
results  are  obtained  when  designing  for  a  maximum  probability 
of  failure  of  0.5  in  each  member;  assuming  a  normal  distribution, 

W   =   2.639  lb. 

Al   =   A3   =   .7883  sq.  in. 

A2   =   .4094  sq.  in. 

The  design  space  is  shown  in  Figure  4.4.   The  lines  of  con- 
stant probability  of  failure  for  SIGH  are  shown.   As  one 
would  expect  the  line  of  constant  0.5  probability  of  failure 
for  SIGH  in  Figure  4.4  is  the  same  as  the  constant  20,000 
psi  line  for  SIGH  in  Figure  4.3.   The  reverse  is  also  true. 
If  the  three  bar  truss  is  optimized  for  a  maximum  probability 
of  failure  of  0.4  the  resulting  safety  factor  for  SIGH  is 
1.03  or  19,361  psi  and  the  minimum  weight  is  2.726  lb.   Using 
19,361  psi  as  the  limit  of  stress  for  a  deterministic  optimi- 
zation yields  approxiamtely  the  same  minimum  weight.   In  both 
cases  the  feasible  design  space  is  the  same  with  the  line 
of  constant  0.4  probability  of  failure  and  constant  19,361 
psi  lying  on  top  of  each  other. 
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Table  I  shows  the  results  of  minimizing  the  weight  of 
the  three  bar  truss  for  various  probabilities  of  failure. 
Table  II  provides  the  deterministic  optimum  designs  for  the 
equivalent  factor  of  safety.   The  difference  in  designs  is 
within  the  numerical  accuracy  of  the  optimization. 

TABLE  I 
THREE  BAR  TRUSS  PROBABILISTIC  DESIGNS 


PROBABILITY 

WEIGHT 

Al    =   A3 

A2 

SIGH 

SAFETY 

OF    FAILURE 

LBS. 

SQ.    IN. 

SQ.    IN. 

PSI 

FACTOR 

.5 

2.6390 

.78827 

.40938 

20000 

1.00 

.4 

2.7260 

.81606 

.41785 

19361 

1.03 

.3 

2.8238 

.84478 

.43438 

18691 

1.07 

.2 

2.9430 

.86887 

.48544 

17940 

1.11 

.1 

3.1208 

.92482 

.50500 

16914 

1.18 

.01 

3.6099 

1.0490 

.64280 

14640 

1.37 

.001 

4.0459 

1.1777 

.71493 

13060 

1.53 

.0001 

4.4770 

1.3212 

.73994 

11792 

1.70 

TABLE  II 
THREE  BAR  TRUSS  DETERMINISTIC  DESIGNS 


SAFETY 

WEIGHT 

Al    =   A3 

A2 

SIGH 

FACTOR 

LBS. 

SQ.    IN. 

SQ.    IN. 

PSI 

1.00 

2.6326 

.77962 

.42752 

20051 

1.03 

2.7212 

.80145 

.45437 

19402 

1.07 

2.8200 

.82773 

.47889 

18726 

1.11 

2.9322 

.87045 

.47014 

18002 

1.18 

3.1180 

.91845 

.52000 

16933 

1.37 

3.5936 

1.0807 

.53693 

14688 

1.53 

4.0258 

1.2033 

.62227 

13110 

1.70 

4.4800 

1.3379 

.69573 

11781 

It  should  be  noted  that  the  correspondence  of  the  proba- 
bility of  failure  and  the  safety  factor  shown  here  is  unique 
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to  this  problem,  in  which  all  C  =  0.1.   In  general  nonlinear 
optimization,  this  one  to  one  relationship  cannot  be  assured. 

C.   EFFECTS  OF  COEFFICIENTS  OF  VARIATION  ON  THE  DESIGN  SPACE 

As  the  coefficients  of  variation  of  the  variables  used  in 
the  design  subroutine  become  smaller,  the  bands  in  the  design 
space  for  the  constrained  values  between  Pf  =  1.0  and  Pf  =  0.0 
becomes  narrower.   Figure  4.5  shows  the  design  space  for  the 
three  bar  truss  in  the  vicinity  of  the  optimum  for  the  case 
where  all  variables  have  coefficients  of  variation,  C  =  0.1. 
Figure  4.6  shows  the  same  space  for  the  case  where  all  coeffi- 
cients of  variation,  C  =  0.001.   In  both  figures  the  lines  of 
constant  objective  and  probabilities  of  failure  for  the 
critical  constrained  variable,  SIGH,  are  shown.   The  line 
of  constant  Pf  =  0.5  is  the  same  for  both  but  the  band  between 
Pf  =  0.9  and  Pf  =  0.1  is  much  narrower  for  the  case  where 
C  =  0.001. 

This  phenomenon  can  cause  numerical  difficulties  in  the 
optimization  process.   Premature  termination  of  optimization 
can  occur  as  the  optimizer  may  be  unable  to  move  down  the 
resulting  'narrow  valley*.   The  absolute  and  relative  termina- 
tion criteria  of  CONMIN  may  have  been  satisfied  but  the  true 
optimum  might  not  have  been  reached. 

If  the  initial  design  is  a  considerable  distance  from 
the  narrow  band  of  changing  probabilities  of  failure,  CONMIN 
may  obtain  zero  gradients  for  the  constraints  which  would 
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Figure  4.5.   THREE  BAR  TRUSS  DESIGN  SPACE  (C  =  0.1) 
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Figure  4.6.   THREE  BAR  TRUSS  DESIGN  SPACE  (C  =  0.001) 
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also  cause  termination  of  the  optimization  process.   A 
routine  to  attempt  to  prevent  this  was  included  in  C0PE19 . 
The  cumulative  distribution  curve  for  the  normal  format  is 
shown  in  Figure  4.7.   The  numerical  integration  routine  used 
to  calculate  the  area  under  the  unit  normal  distribution 
curve  yields  a  probability  of  failure,  Pf  of  1.0000000  for 
a  6  of  -5.2999821  and  a  Pf  of  0.0000000  for  a  6  of  5.5000086 
It  was  modified  to  "widen"  the  zone  between  Pf  =  1.0  and 
Pf  =  0 .   A  polynomial  (a  line  equation)  was  placed  at  each 
end  of  the  cumulative  distribution  routine  to  yield  a  Pf 
of  1.0  at  3  equals  -25.0,  decreasing  to  a  Pf  of  0.9999866 
at  -4.2  and  a  Pf  of  0.00001335144  at  3  equals  4.2  decreasing 
to  a  Pf  of  0.0  at  25.0.   This  is  shown  in  Figure  4.8.   This 
provides  the  same  optimization  results  within  the  accuracy 
of  the  program,  but  avoids  the  problem  of  zero  gradients  of 
the  constrained  variables. 

An  additional  approximation  was  used  to  determine  the 
probability  of  failure  when  the  coefficients  of  variation 
concerned  are  zero  or  nearly  so.   The  probability  that  X  is 
greater  than  Y  is,  in  the  standard  format 


■B[-B]   =   V"  Sf1  (48) 


where 


yz   =   ^y  "  Ux 


(49) 
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Figure  4.8.   MODIFIED  CUMULATIVE  NORMAL  DISTRIBUTION  CURVE 
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and 


°Z        =    J°l    +  ay  (50) 


Mathematically,  if  a   =  0 ,  then  one  of  three  situations  can 
occur. 


1.  X  =  Y,  8  =  0,  then  Pf  =  0.5; 

2.  X  >  Y,  8  ■*    -°°,  then  Pf  =  1.0; 

3.  X  <  Y,  8  ■+  °°,  then  Pf  =  0.0. 

This  is  shown  graphically  in  Figure  4.9.   In  order  to  provide 
CONMIN  with  a  smooth  function  with  which  to  work,  a  polynomial 
was  used  to  connect  the  Pf  =  1.0  and  Pf  =  0.0  in  Figure  4.9. 
This  is  shown  in  Figure  4.10.   The  interpolating  polynomial 
is  used  when  8/  the  modified  safety  index,  is  within  one 
unit  of  0.   Where  8  is  defined  to  be 


8   =     Y  "  X  (51) 

v/YTT7r 


A  similar  approximation  was  developed  for  the  lognormal 
mo  de  1 . 

Experience  has  shown  that  these  modifications  dramatically 
improve  the  numerical  stability  of  the  optimization  process 
without  significantly  affecting  the  mathematical  result. 
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V.   COMPUTER  PROGRAM 

A.  COPES  INPUT  DATA 

Standard  deterministic  optimization  requires  that  COPES 
data  blocks  A-0  and  V  be  used.   For  probabilistic  optimiza- 
tion the  same  data  is  entered  in  these  blocks  with  the 
exception  of  NCALC  and  IPROB  in  data  block  B.   Here  an  NCALC 
=  7  would  be  entered  to  indicate  probabilistic  optimization 
was  to  be  performed  and  an  IPROB  =  1  or  2  to  indicate  which 
probability  model  is  to  be  used.   Additional  data  blocks  are 
required  for  probabilistic  optimization.   Appendix  A  con- 
tains detailed  input  instructions  for  these  data  blocks  as 
well  as  an  example  problem  including  sample  input  data. 

B.  OPTIMIZATION  DATA  FLOW 

The  interaction  between  COPES,  CONMIN ,  and  the  analysis 
subroutine  will  be  the  same  for  probabilistic  optimization 
as  it  is  in  standard  deterministic  optimization  with  one 
exception.   In  the  optimization  process,  each  time  COPES 
requires  the  determination  of  a  constraint  vector  it  will 
call  COPE  19.   C0PE19  will  call  the  analysis  subroutine  as 
many  times  as  necessary  to  determine  the  required  partial 
derivatives  by  finite  difference  steps.   Once  these  are 
obtained,  the  constraint  vector  is  determined  and  provided 
to  COPES.   A  flow  diagram  for  probabilistic  optimization 
is  provided  in  Figure  5.1. 
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C.   COPE 19  SUBROUTINE 

The  main  computer  subroutine  developed  in  this  investi- 
gation is  C0PE19.   Three  additional  subroutines,  COPE20 
through  COPE22,  were  also  developed  to  perform  specific 
tasks  required  by  C0PE19 .   The  operation  and  flow  of  these 
subroutines  are  discussed  below. 

1.   Input 

The  following  is  provided  to  C0PE19: 

a.  The  COPES  control  parameter  ICALC. 

b.  Which  probability  model  to  use,  either  normal 
or  lognormal. 

c.  The  number  of  probability  variables  which  make 
up  the  constraint  equations  (IVAR) . 

d.  The  global  locations  in  the  ANALIZ  subroutine 
common  block  of  each  IVAR. 

e.  Coefficient  of  variation  of  each  IVAR,  assumed 
constant  throughout  the  optimization. 

f.  The  number  of  constrained  variables. 

g.  The  global  location  in  the  ANALIZ  subroutine 
common  block  of  each  of  the  constrained 
variables . 

h.   The  upper  and  lower  limits  imposed  on  each  of 

the  constrained  variables, 
i.   The  allowed  probability  of  failure  at  the  upper 

and  lower  limits  of  each  constrained  variable, 
j .   The  coefficient  of  variation  of  the  constrained 

variable  limits. 
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k.   The  variable  MGPAD  which  determines  how  often 

the  constrained  variable  gradients  are  calculated 
2.   Determination  of  Constraint  Vector 

When  COPES  requires  a  constraint  vector  for  a  given 
design  it  will  call  C0PE19  with  ICALC  =  2.   C0PE19  will  then 
calculate  the  partial  derivatives  of  each  constrained  varia- 
ble with  respect  to  the  variables  which  make  up  the  con- 
straint equations.   The  subroutine  COPE20  controls  the  finite 
difference  calculations.   In  the  example  in  Section  IV. A, 
the  cantilevered  beam,  it  would  calculate: 


3SB      6P 


3L 

BH2 

3SB 

6L 

9P 

BH2 

3SB 

12PL 

3H 

BH3 

3SB 

6PL 

3B 

2  2 
B  H 

(52) 


(53) 


(54) 


(55! 


C0PE19  will  then  calculate  the  standard  deviation  (Equation 
44)  of  the  constrained  variable,  in  this  case  SB,  as  dis- 
cussed in  Section  III.B  and  Section  IV. A.   Now  the  safety- 
index  is  calculated  and  subroutine  C0PE21  called  to  deter- 
mine the  probability  of  failure.   When  required,  C0PE21  will 
call  subroutine  COPE22  to  perform  numerical  integration  of 
the  normal  density  function.   The  constraint  value  to  be 
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stored  in  the  G  vector  is  then  determined  (Equation  6) . 
The  above  would  be  performed  for  each  constrained  variable 
and  the  resulting  G  vector  provided  to  COPES. 

D.   REDUCTION  OF  FUNCTION  EVALUATIONS 

In  order  to  determine  the  gradients  of  the  constrained 
variables  finite  difference  steps  are  taken  resulting  in  a 
significant  increase  in  the  number  of  analyses  performed  in 
the  optimization  process.   In  many  problems  these  gradients 
remain  essentially  the  same  throughout  the  optimization. 
Therefore,  it  is  sometimes  possible  to  reduce  the  frequency 
with  which  the  gradients  are  calculated  and  achieve  approxi- 
mately the  same  result.   The  input  variable  code  LGRAD  is 
included  as  an  option  in  probabilistic  optimization  to  allow 
this  choice.   LGRAD  has  the  following  meanings: 

0:   Calculate  gradients  each  time  C0PE19  is  called. 

1:   Calculate  gradients  at  the  beginning  of  each 
CONMIN  iteration. 

2;   Calculate  gradients  only  at  the  beginning  of 
optimization. 

Table  III  demonstrates  the  result  of  various  LGRAD  selections 
when  designing  the  three  bar  truss  for  an  allowed  probability 
of  failure  of  0.50,  where  all  coefficients  of  variations, 
C  =  0.1. 

Probabilistic  optimization  was  performed  using  the  three 
LGRAD  values  on  numerous  additional  test  cases.   In  almost 
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TABLE  III 
REDUCTION  OF  FUNCTION  EVALUATIONS:   THREE  BAR  TRUSS  DESIGN 


LGRAD    WEIGHT    Al  =  A3      A2       NUMBER  OF 

LBS.     SQ.  IN.    SQ.  IN.    FUNCT  EVALS 


Deterministic 

2.6326 

.77962 

.42752 

35 

Probabilistic 

0 

2.6390 

.78827 

.40938 

343 

Prcbabii is tic 

1 

2.6389 

.77386 

.45011 

244 

Probabilistic 

2 

2.6379 

.78972 

.40427 

56 

every  case  the  optimization  results  were  essentially  the  same 
for  LGRAD  =  0  and  1.   The  composite  driveshaft  design  in 
Chapter  V  provides  an  excellent  example  showing  the  function 
evaluation  reductions  when  performing  a  relatively  complex 
design.   Table  IV  demonstrates  the  results  of  reducing  the 
frequency  of  constrained  variable  gradient  calculations 
when  designing  the  composite  driveshaft  for  minimum  weight, 
using  the  normal  distribution  model.   For  this  particular 
problem,  calculation  of  constrained  variable  gradients  only 
at  the  beginning  of  optimization  was  insufficient  to  obtain 
results  near  the  optimum. 

TABLE  IV 
REDUCTION  OF  FUNCTION  EVALUATIONS:   DRIVESHAFT  DESIGN 


PAZ 

WEIGHT 
L3S. 

0 
1 
2 

9.0203 
9.0205 
27.216 

NUMBER  OF  FUNCTION 
EVALUATIONS 

4011 

1344 

89 
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VI.   NUMERICAL  EXAMPLES 

The  following  examples  are  presented  to  demonstrate  the 
capabilities  of  the  computer  subroutines  developed  in  this 
investigation.   In  both  cases  an  existing  ANALIZ  subroutine 
was  used  after  minor  modifications.   For  these  two  shaft 
designs  the  subroutine  from  [Ref .  6] ,  was  used.   The  ANALIZ 
subroutine  was  altered  to  include  the  required  probability 
variables  in  the  global  common  block  and  to  remove  the 
factor  of  safety  calculations  which  were  originally  used  as 
constraints . 

In  each  example  a  similar  deterministic  design  was  per- 
formed to  give  the  reader  a  better  feel  for  the  results. 
The  deterministic  designs  are  not  intended  to  duplicate  the 
probabilistic  design  results. 

The  mean  values  used  in  the  designs  are  from  [Ref.  6]. 
The  coefficients  of  variation  are  reasonable  assumed  values 
used  for  demonstration  purposes. 

A.   STEEL  DRIVESHAFT 

Design,  for  minimum  weight,  a  steel  driveshaft  to  trans- 
mit 150.0  horsepower  at  300.0  RPM.   The  loading  is  presented 
in  Figure  6.1.   The  shaft  is  designed  against  failure  in 
strain,  in  torsional  buckling,  vibration  frequency  and  a 
maximum  deflection  of  0.0  5  inches. 

The  design  variables  are  the  shaft  thickness  and  the 
shaft  inside  diameter.   The  thickness  is  to  be  maintained 
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Figure  6.1.   DRIVESHAFT  LOADING 
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between  0.01  and  2.0  inches  and  the  inside  diameter  between 
0.5  inches  and  5.0  inches.   The  initial  design  is  an  inside 
diameter  of  2.0  inches  and  a  thickness  of  1.0  inch.   The 
remaining  data  for  the  analysis  is  presented  in  Table  V. 

TABLE  V 
INPUT  DATA  FOR  STEEL  DRIVESHAFT 


HORSEPOWER 
SPEED 

AXIAL  LOAD 
LENGTH 

YOUNGS  MODULUS 
SHEAR  MODULUS 
POISSONS  RATIO 
SPECIFIC  WEIGHT 
STRAIN  LIMIT 
MAX  DEFLECTION 
THICKNESS 
INSIDE  DIAMETER 
FAILURE  DUE  TO: 

TORSIONAL  BUCKLING 

FREQUENCY 


MEAN 

COEFFICIENT 

VALUE 

OF 

VARIATION 

150. 

HP 

.08 

300. 

RPM 

.08 

50. 

LB. 

.18 

120. 

IN. 

.015 

30,000,000 

PSI 

.038 

11,538,000 

PSI 

.042 

.3 

.026 

.2  82 

LB/CU.  IN. 

.01 

.001 

IN. /IN. 

.01 

.05 

IN. 

.05 

.015 

.015 

.21 
.06 

Allowing  for  maximum  probabilities  of  failure  of  0.01  in 
strain,  torsional  buckling,  deflection  and  vibration  fre- 
quency; a  minimum  weight  of  54.15  lbs  is  achieved  for  the 
normal  model  and  52.88  lbs  for  the  lognormal  model.   The 
initial  design  summary  is  presented  in  Table  VI  and  the 
final  designs  summaries  and  COPES  optimization  results  in 
Tables  VII  through  XII. 

Designing  for  a  factor  of  safety  against  failure  of  2.0 
for  strain,  torsional  buckling,  vibration  frequency,  and 
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TABLE    VI 
STEEL    DHIVESHAFT    INITIAL    DESIGN    SUMMARY 


STEEL    DRIVESHAFT    OUTPUT 


NUMBER    OF    PLYS  =  1 

NUMBER    OF    MATERIAL    TYPES       =  1 

NUMBER    OF    LOAD    CONDITIONS    =  1 

ECCENTRICITY  =  0.0 


DIMENSIONS 

PLY  THICKNESS         PERCENT  DIAMETER         THETA 

INSIDE    0.20000E+01 
1    0. 10G00E+01       100-00         0.40000E+01         0.0 


STIFFNESS 
AE    =    0.31071E+09 
EI    =    0.38838E+09 
GJ    =    0.27187E+09 


LOADS: 

L~  T  M  F  H.  °  RPM 

l""0.31513E+05       O'.O  0.50000E  +  02    0.15000E+03    0.30000E+03 


LOAD    CONDITION  1 

PLY             PPL                  S.F.  ECT                  S.F.                   EPLT                  S.F. 

1    0.16092E-06    100!o5  oTo                 100.00       0.23182E-03         5*.  61 

CRITICAL    SPEED  =  0.15543E  +  0'4 

MAXIMUM    DEFLECTION    =  0.19195E-01 

WEIGHT  =  0  .  318  9  3E+0  3 

VOLUME  =  0  .  11310E  +  0U 
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TABLE   VII 
STEEL    DRIVESHAFT    COPES    STANDARD    OUTPOT:    NORMAL 


OPTIMIZATION    RESULTS 

OBJECTIVE    FUNCTION 

GLOEAL    LOCATION  3  FUNCTION    VALUE    0.54154E+02 

DESIGN    VARIABLES 

D.    V.  GLOBAL  LOWER  UPPER 

ID  NO.  VAR.     NO.  BOUND  VALUE  BOUND 

1  1  11  0.10000E-01  0.99883E-01  0.20000E+01 

2  2  2  0.50000E+00  0.50000E+01  0.50000E+01 


DESIGN   CONSTRAINTS 

DETERMINISTIC 

GLOBAL  LOWER  MEAN  UPPER 

ID  VAR.    NO.  BOUND  VALUE  BOUND 

1  51  -0.1GOO0E-02  0.94774E-06  0.10000E-02 

3  52  -0.10000E-02  0.0  0.10000E-02 

5  53  -0.10000E-02  0.68213E-03  0.10000E-02 

7  54  0.31513E*05  0.243'40E  +  06  0.11000E+16 

8  56  0.30000E+03  0.25068E+04  0.11000E+16 

9  57  -0.11000E+16  0.72076E-02  0.50000E-01 
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TABLE   VIII 
STEEL    DRIVESHAFT    COPES    PROBABILISTIC    OUTPUT:     NORMAL 

PROBABILISTIC 


MEAN 

STANDARD 

COEFFICIENT 

ID 

VALUE 

DEVIATION 

OF    VARIATION 

1 

0.947742-06 

0. 17535-06 

0.  1753E-01 

3 

0.0 

0.0 

0.0 

5 

0.68213E-03 

0.6535E-04 

0. 9581E-01 

7 

0.24340E+06 

0.1364E+05 

0. 5606E-01 

8 

0.25068S+04 

0.9634E+02 

0. 38U3E-01 

9 

0.72076E-02 

0.5707E-03 

0.7918E-01 

PROBABILITY    OF    FAILURE  PROBABILITY    OF    FAILURE 

AT    LOWER    BOUND  AT    UPPER    BOUND 


ID 

ALLOWED 

CALCULATED 

ALLOWED 

CALCULATED 

1 

0.  1000E-01 

0.0 

0.  1000S-01 

0.0 

3 

0.  1000E-01 

0.0 

0.  1000E-01 

0.0 

5 

0.  1000E-01 

0.0 

0.  1000E-01 

0. 1000E-01 

7 

0.  1000E-01 

0.0 

0. 1000E+01 

0.0 

8 

0.1000E-01 

0.0 

0.  1000E  +  01 

0.0 

9 

0.1000E+01 

0.0 

0.  1000E-01 

0.0 
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T1BLE    IX 
STEEL    DRIVESHAFT    DESIGN    SUMMARY:    NORMAL 


STEEL    DRIVESHAFT   OUTPUT 


NUMBER    OF    PLYS  =  1 

NUMBER    OF    MATERIAL    TYPES      =  1 

NUMBER    OF    LOAD   CONDITIONS    =  1 

ECCENTRICITY  =  0.0 

DIMENSIONS 

PLY  THICKNESS         PERCENT         DIAMETER         THEIA 

INSIDE    0.50000E+01 
1    0.99883S-01       100.00         0.51993E+01         0.0 


STIFFNESS 
AE    =    0.52757E+08 
EI    =    0.17158E+09 
GJ    =    0.12011E+09 


LOADS: 

L.C.  T  M  F  H.P.  RPM 

1    0.31513E+05      0.0  0.53000E+02      0.15000E+03       0.300002+03 


LOAD    CONDITION  1 

FLY            E?L                 S.F.  EPT                 S.F.               EPLT                    S.F. 

1      0.94774E-06    100.00      0.0  100.00       0.68213E-03         1.91 

CRITICAL    SPEED  =    0.25068E  +  0!* 
MAXIMUM    DEFLECTION    =    0.72076E-02 

WEIGHT  =    0  .54154E+02 

VOLUME  =    0  .192  04E+03 
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TABLE    X 


STEEL    DRIVESHAFT    COPES    STANDARD    OUTPUT:    LOGNORMAL 


OPTIMIZATION    RESULTS 


OBJECTIVE    FUNCTION 
GLOBAL    LOCATION  3 


FUNCTION    VALUE    0.52882E+02 


DESIGN    VARIABLES 


ID 

1 

2 


D.     V. 
NO. 
1 
2 


GLOBAL 
VAR.    NO. 
11 
2 


LOWER 
BOUND 
0.10000E-01 
0.5030 OE+00 


VALUE 
0. 97581E-01 
0.50000E+01 


UPPER 
BOUND 
0.20Q00E+01 
0.50000E+01 


DESIGN    CONSTRAINTS 


DETERMINISTIC 


GLOBAL 

LOWER 

MEAN 

UPPER 

ID 

VAR. 

NO. 

BOUND 

VALUE 

BOUND 

1 

51 

-0.10000E-02 

0.97054E-06 

0. 

10000E-02 

3 

52 

-0.10000E-02 

0.0 

0. 

,  10000F-02 

5 

53 

-0.10000S-02 

0.59856E-03 

0. 

10000E-02 

7 

54 

0.31513E+05 

0.22961E+06 

0. 

11000E+16 

8 

56 

0.30000E+03 

0.25056E+0U 

0. 

11000E+16 

9 

57 

-0.11000  E+16 

0.72144E-02 

0, 

.50C00E-01 
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TABLE    XI 
STEEL    DRIVESHAFT    COPES    PROBABILISTIC    OUTPUT:     LOGNORMAL 

PROBABILISTIC 


MEAN 

STANDARD 

COEFFICIENT 

ID 

VALUE 

DEVIATION 

OF    VARIATION 

1 

0.97054E-06 

0. 1796S-06 

0.  1796E-01 

3 

0.0 

0.0 

0.0 

5 

0. 69856E-03 

0.6692E-04 

0.9580E-01 

7 

0.22961E+06 

0.  1283E+05 

0.5611E-01 

8 

0.  25056E+04 

0.9630E+02 

0.38432-01 

9 

0.72144E-02 

0.5713E-03 

0.7918E-01 

PROBABILITY    OF    FAILURE  PROBABILITY    OF    FAILURE 

AT    LOWER    BOUND  AT    UPPER     BOUND 


ID 

ALLOWED 

CALCULATED 

ALLOWED 

CALCULATED 

1 

0. 1000E-01 

0.0 

0.  10003-01 

0.0 

3 

0.  1000E-01 

0.0 

0.  1000E-01 

0.0 

5 

0. 1000E-01 

0.0 

0. 1000E-01 

0.9993E-02 

7 

0.  1000E-01 

0.0 

0.  1000E+01 

0.0 

8 

0. 1000E-01 

0.0 

0.  1000S  +  01 

0.0 

9 

0.  1000E+01 

0.0 

0. 1000E-01 

0.0 
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TABLE   XII 
STEEL    DRIVESHAFT    DESIGN    SUMMARY:     LOGNORMAL 


STEEL    DRIVESHAFT   OUTPUT 


NUMBER    OF    PLYS  =  1 

NUMBER    OF    MATERIAL    TYPES  =  1 

NUMBER    OF    LOAD    CONDITIONS  =  1 

ECCENTRICITY  =  0.0 


DIMENSIONS 

PLY  THICKNESS         PERCENT         DIAMETER  THSTA 

INSIDE    0.50000E+01 
1    0.97581E-01       100.00         0.51952E+01         0.0 


STIFFNESS 
AE    =    0.51518E+08 
EI    =    0. 16740E+09 
GJ    =    0.1 1718E+09 


LOADS : 

L.C.  T  M  F  H.P.  RPM 

1       0.31513E+05      0.0        0.50000E+02      0.15000E+03       0.30000E+03 


LOAD    CONDITION  1 

PLY  EPL  S.F.  EPT  S.  F.  EPLT  S.F. 

1       0.97054E-06    100.00  0.0  100.00         0.69856E-03         1.86 

CRITICAL    SPEED  =    0.25056E  +  0'4 

MAXIMUM    DEFLECTION  =    0.72144E-02 

WEIGHT  =    0 .52882E+02 

VOLUME  =    0 .  18752E  +  03 
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deflection  an  optimum  of  73.78  lbs  is  obtained.   The  COPES 
optimization  results  and  final  design  summary  are  contained 
in  Tables  XIII  and  XIV. 

B.   COMPOSITE  DRIVESHAFT 

Design,  for  minimum  weight  a  120  in.  long,  four  ply 
graphite  epoxy  driveshaft  to  transmit  150.0  horsepower  at 
300  RPM.   Figure  6.1  shows  the  shaft  loading  condition. 
The  shaft  is  designed  against  failure  in  transverse,  logitu- 
dinal,  and  shear  strain;  torsional  buckling,  vibration  fre- 
quency and  a  maximum  deflection  of  0.0  5  inches.   Table  XV 
contains  the  initial  design. 

The  design  variables  are  the  thickness  of  each  ply,  the 
orientation  of  each  ply  and  the  inside  diameter.   Addition- 
ally, the  thickness  of  ply  2  and  ply  3  must  remain  equal  to 
each  other  and  the  orientation  of  ply  3  equals  the  negative 
of  ply  2.   The  following  additional  constraints  are  placed 
on  the  design  variables:   The  inside  diameter  must  remain 
between  0.5  inches  and  5.0  inches  and  the  thickness  of  any 
given  ply  between  0.01  and  0.5  inches.   Ply  2  is  allowed  to 
vary  between  0.0  and  9  0.0  degrees.   The  orientation  of  plys 
1  and  4  remain  constant.   The  remaining  input  data  is  pre- 
sented in  Table  XVI. 

Allowing  for  maximum  probabilities  of  failure  of  0.01  in 
strain,  torsional  buckling,  deflection,  and  vibration  fre- 
quency; an  optimum  weight  of  9.02  lbs  is  achieved  using  the 
normal  distribution  model  and  10.28  lbs  for  the  lognormal 
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TABLE   XIII 
STEEL    DHI7ESHAPT    COPES    OUTPUT :    DETERMINISTIC 


OPTIMIZATION    RESULTS 

OBJECTIVE    FUNCTION 

GLOBAL    LOCATION  3  FUNCTION    VALUE    0.73783E+02 

DESIGN    VARIABLES 


D.    V. 
ID         NO. 

1  1 

2  2 

GLOBAL 
VAR.    NO. 
11 
2 

LOWER 

BOUND 
0.10000E-01 
0.50000E+00 

VALUE 
0. 13515E+00 
0. 50000E+01 

UPPER 

BOUND 
0.20000E+01 
0.50000E+01 

DESIGN    CONSTRAINTS 

DETERMINISTIC 

ID 

1 

3 

5 

7 

8 

9 

GLOBAL 
VAR.    NO. 

51 

52 

53 

54 

56 

57 

LOWER 

BOUND 

-0.50000E-03 

-0.50000E-03 

-0.50000E-03 

0 .63025  E  +  05 

0.60000E+03 

-0. 11 000 E* 16 

MEAN 

VALUE 
0.59561E-06 
0.0 

0.50035E-03 
0.51838E+06 
0.25247E+04 
0.71044E-02 

UPPER 

BOUND 
0.50000E-03 
0.50000F-03 
0.50000E-03 
0.  11000E+16 
0.  11000E+16 
0.20000E-01 
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TABLE   XIV 
STEEL    DRIVESHAFT    DESIGN    SUMMARY:     DETERMINISTIC 

STEEL    DRIVESHAFT   OUTPUT 

NUMBER    OF    PLYS  =  1 

NUMBER    OF    MATERIAL    TYPES  =  1 

NUMBER    OF    LOAD    CONDITIONS  =  1 

ECCENTRICITY  =  0.3 

DIMENSIONS 

PLY  THICKNESS         PERCENT  DIAMETER         THETA 

INSIDE    3.50000E  +  01 
1    0.  13515E  +  00       100.00         0.52703E+01         0.0 

STIFFNESS 
AE    =    0.71879S+08 
EI    =    0.23709E+09 
GJ    =    0.16597E+09 


LOADS: 

L.C.  T  M  F  H.P.  RPM 

1       0. 315132+05      0.0         0.50000E+02      3.150Q0E+03       0.300C0E+03 


LOAD    CONDITION  1 

PLY            EPL                 S    F  pdj                 S    p                    SPLT                    S.F. 

1       0.69561E-06    100.00  0.0                 100lo6         0.500352-03         2*.60 

CRITICAL    SPEED  =    0  .252 4 72  +  0'4 
MAXIMUM    DEFLECTION    =    3.710U4E-02 

WEIGHT  =    0  .737832+02 

VOLUME  =    0  .26164E+03 
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TABLE  XV 
COMPOSITE  DRIVESHAFT  INITIAL  DESIGN 


PLY 
PLY 
PLY 
PLY 


1 
2 

3 
4 


THICKNESS 

ORIENTATION 

IN. 

DEGREES 

.25 

0.0 

.25 

20.0 

.25 

-20.0 

.25 

90.0 

INSIDE  DIAMETER    2.0  IN. 

TABLE  XVI 
INPUT  DATA  FOR  COMPOSITE  DRIVESHAFT 


HORSEPOWER 

SPEED 

LENGTH 

AXIAL  LOAD 

LONGITUDINAL  MODULUS 

TRANSVERSE  MODULUS 

SHEAR  MODULUS 

MAJOR  POSSIONS  RATIO 

STRAIN  LIMITS: 

LONGITUDINAL  COMPRESSIVE 
LONGITUDINAL  TENSILE 
TRANSVERSE  COMPRESSIVE 
TRANSVERSE  TENSILE 
MAX  SHEAR 

SPECIFIC  WEIGHT 

MAX  DEFLECTION 

PLY  ORIENTATION 

PLY  THICKNESS 

INSIDE  DIAMETER 

FAILURE  DUE  TO: 

TORSIONAL  BUCKLING 
FREQUENCY 


MEAN 

COEFFICIENT 

VALUE 

OF 

VARIATION 

150. 

HP 

.08 

300. 

RPM 

.08 

120. 

IN. 

.015 

50. 

LB. 

.18 

21,000,000 

PSI 

.14 

1,700,000 

PSI 

.14 

650,000 

PSI 

.14 

.21 

.08 

-.00857 

IN. /IN. 

.15 

.00857 

IN. /IN. 

.15 

-.0176 

IN. /IN. 

.15 

.00471 

IN. /IN. 

.15 

.0184 

IN. /IN. 

.15 

.056 

LB/CU. 

IN. 

.05 

.05 

IN. 

.05 

.05 

.015 

.015 

.20 
.06 
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model.   The  ANALIZ  subroutine  design  summary  for  the  initial 
design  is  presented  in  Table  XVII.   Tables  XVIII  through 

XXIII  contain  the  design  summaries  and  final  COPES  optimiza- 
tion results. 

Designing  for  factors  of  safety  against  failure  of  2.0 
for  strain,  torsional  buckling,  deflection,  and  vibration 
frequency  an  optimum  weight  of  10.23  lbs  is  obtained.   Tables 

XXIV  and  XXV  present  the  deterministic  design  results. 
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TABLE   XVII 
COMPOSITE    DRIVESHAFT    INITIAL    DESIGN    SUMMARY 

COMPOSITE    DRIVESHAFT    OUTPUT 

NUMBER    OF    PLYS  =  4 

NUMBER    OF    MATERIAL    TYPES  =  1 

NUMBER    OF    LOAD   CONDITIONS  =  1 

ECCENTRICITY  =  3.0 

DIMENSIONS 

PLY         THICKNESS  PERCENT         DIAMETER  THETA 
INSIDE    0.20000E+01 

1  0.25000E+00  25.00         0. 25000E+0  1  0.0 

2  0.25000E+00  25.00         0.30003E+01  0.20000E+02 

3  0.25000E+00  25.00         0.35000E+01  -0.20000E+02 

4  0. 25000E+00  25.00         0.40003E+01  0 . 930 00E+02 

STIFFNESS 
AE    =    0. 12143E+09 
EI    =    0.12426E+09 
GJ    =    0.37225E+08 


LOADS: 

L   C                  T                          M  F                               HP                             RPM 

1    "0-31513S+05      0.0  0.50000E+02      0.15000S+03      O.3O0OOE+O3 

LOAD    CONDITION  1 

PLY            E^L                    S    F  EPT                 S.p.                  EPLT               S.17. 

1  0.41176E-06    1OOl0O  0.0                         IOOI60       0.10582E-02    17~.39 

2  0.40847E-03       20.98  -0.40806E-03      U  3. 1 3       0.97246E-03    18.92 

3  -0. 47576E-03       18.01  0.47617E-03         9.89       0.11351E-02    16.21 

4  0.53159E-09    100.00  3.41123E-06    100.00    -0.16931E-02    10.87 

CRITICAL    SPEED  =    0.19725E+04 
MAXIMUM    DEFLECTION    =    0.11746E-01 

WEIGHT  =    0 .63334E+02 

VOLUME  =    0  .11310E  +  04 
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TABLE  XVIII 
COMPOSITE  DRIVESHAFT  COPES  STANDARD  OUTPUT:  NORMAL 

OPTIMIZATION  RESULTS 

OBJECTIVE  FUNCTION 

GLOBAL  LOCATION     3      FUNCTION  VALOE  0. 90203E+01 


DESIGN 

VARIABLES 

D. 

V. 

GLOBAL 

LOWER 

UPPER 

ID         NC 

i. 

VAR.     NO. 

BOUND 

VALUE 

BOUND 

1             1 

11 

0.10000E-01 

0. 

10000E-01 

0.50000E+00 

2           2 

12 

0.10000E-01 

3. 

10212E+00 

0.50000E+00 

3           2 

13 

0.10000E-01 

0. 

10212E+00 

0.50000E+00 

4           3 

14 

0.10000E-01 

0. 

10000E-01 

0.50000E+09 

5            4 

32 

0.0 

0. 

34968E+02 

0.90000E+02 

6            4 

33 

0.0 

-0. 

34968E+02 

-0.90000E+02 

7            5 

2 

0.50000E+00 

0. 

16311E+01 

0.50000E+01 

DESIGN 

CONSTRAINTS 

DETERMI NISTIC 


GLOBAL 

LOWER 

MEAN 

UPPER 

ID 

VAR.    NO. 

BOUND 

VALUE 

BOUND 

1 

51 

-0.  85700  E-02 

0.35785E-05 

0. 

.85700E-02 

3 

52 

-0.17600E-01 

0.0 

0. 

■47100E-02 

5 

53 

-0.18400 E-01 

0.47679E-02 

0. 

,  18400E-01 

7 

54 

-0.85700E-02 

0.25105E-02 

0. 

.85700E-02 

9 

55 

-0.17600E-01 

-0.25069E-02 

0. 

,47100E-02 

11 

56 

-0.18400E-01 

0. 18288E-02 

0. 

18400E-01 

13 

57 

-0.85700E-02 

-0.27745E-02 

0. 

■85730E-02 

15 

58 

-0.17630E-01 

0.27781E-02 

0. 

47100E-02 

17 

59 

-0.18400E-01 

0.20319E-02 

0. 

,  18400E-01 

19 

60 

-0.35700E-02 

0. 18741E-08 

0. 

85730E-02 

21 

61 

-0.17600E-01 

0.35766E-05 

0. 

,47100E-02 

23 

62 

-0.18403E-01 

-0.59688E-02 

0. 

,  18400E-01 

25 

63 

0.31513E+05 

0.55160E+05 

0. 

,  11000E+16 

26 

65 

0.30000E+03 

0. 1  1715E  +  04 

0, 

,  11000E+16 

27 

66 

-0.1100  0  E+ 16 

0.34815E-01 

0. 

.50000E-01 

64 


TABLE   XIX 
COMPOSITE    DRIVESHAFT   COPES    PROBABILISTIC    OUTPUT:    NORMAL 

PROBABILISTIC 

MEAN 

ID  VALUE 

1  0.35785E-05 

3  0.0 

5  0.47679E-02 

7  0.25105E-02 

9  -0.25069E-02 

11  0.18288E-02 

13  -0.27745E-02 

15  Q.27781E-02 

17  0.20319E-02 

19  0.18741E-08 

21  0.35766E-05 

23  -0.59688E-02 

25  0.55160E+05 

26  0.11715E+04 

27  9.34815E-01 

PROBABILITY    OF    FAILURE  PROBABILITY    OF    FAILURE 

AT    LOWER    BOUND  AT    UPPER    BOUND 


STANDARD 

COEFFICIENT 

DEVIATION 

OF    VARIATION 

0.4929E-06 

0.49292-01 

0.0 

0.0 

0.7317E-03 

0. 1535E+00 

0.3838E-03 

0. 1529E+00 

0.38342-03 

0. 15302+CO 

0. 4363E-03 

0.23862+00 

0.4221E-03 

0. 1521E+00 

0.4225E-03 

0. 1521E+00 

0.4938E-03 

0. 2430E+00 

0.U687E-03 

0. 10002+02 

0.4687E-93 

0. 10002+02 

0.9181E-03 

0. 1538E+00 

0.7849E+04 

0. 1423E+00 

0.9463E+02 

0.8078E-01 

0.6032E-02 

0. 1733E+00 

ID 

ALLOWED 

CALCULATED 

ALLOWED 

CALCULATED 

1 

0.  1000E-01 

0.0 

0.  1000E-01 

0.0 

3 

0.  1000E-01 

0.0 

0.  1000E-01 

0.0 

5 

0.  1000E-01 

0.0 

0.  1000E-01 

0. 9537E-06 

7 

0.1000E-01 

0.0 

0. 1000E-01 

0.3338E-05 

9 

0.  1000E-01 

0.0 

0.  1000E-01 

0.0 

11 

0.  1000E-01 

0.0 

0.  1000E-01 

0.0 

13 

0.  1000E-01 

0.9537E-05 

0.  1000E-01 

0.0 

15 

0.  10002-01 

0.0 

0.  1000E-01 

0. 9468E-02 

17 

0.  1000E-01 

0.0 

0. 1000E-01 

0.0 

19 

0.  1000E-01 

0.0 

0.  1000E-01 

0.0 

21 

0.  1000E-01 

0.0 

0.  1000E-01 

0.0 

23 

0.  1000E-01 

0.1001E-04 

0. 1000E-01 

0.0 

25 

0.  1000E-01 

0.9410S-02 

0. 1000E+01 

0.0 

26 

0.  1000E-01 

0.0 

0. 1000E+01 

0.0 

27 

0. 1000E+01 

0.0 

0. 1000E-01 

0.  1002E-01 
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TABLE    XX 
COMPOSITE    DRIVESHAFT    DESIGN    SUMMARY:    NORMAL 

COMPOSITE    DRIVESHAFT    OUTPUT 

NUMBER    OF    PLYS  =  4 

NUMBER    OF    MATERIAL    TYPES  =  1 

NUMBER    OF    LOAD    CONDITIONS  =  1 

ECCENTRICITY  =  0.0 

DIMENSIONS 

PLY         THICKNESS  PERCENT         DIAMETER  THEIA 
INSIDE    0.15  81 1E  +  01 

1  0.10000E-01  4.46         0.17011E+01  0.0 

2  0.  10212S+00  45.54         0.19054E+01  0.34968E+02 

3  0.  10212E+00  45.54          0.21095E+01  -0.34968E+02 

4  0.  10000E-01  4.45         0.21295E+01  0.90000E+02 

STIFFNESS 
AE    =    0.  13972E+08 
EI    =    0.631 16E+07 
GJ    =    0.56217E+07 


LOADS  : 
1    "0.31513E+05         0.0      0.500000E  +  02    0.15000S+03    0.30000E+03 


LOAD 

CONDITION            1 

PLY 

EPL                 S.F.                  SPT                 S. 

,F. 

EPLT 

S. 

,F. 

1       0, 

■35785E-05    100.00      0.0                         100. 

00 

0. 

.47579E- 

■02 

3. 

.86 

2      0. 

■  25105E-02         3.41    -0.25069E-02         7. 

02 

0. 

,  18288E- 

■02 

10, 

.06 

3   -0. 

,  27745E-02         3.09      0.27781E-02         1. 

70 

0. 

.20319E- 

■02 

9. 

.06 

4      0, 

,  18741E-08    100.00      0.35766E-05    100. 

CRITICAL    SPEED              =    0.11715E+04 
MAXIMUM    DEFLECTION    =    0.34815E-01 
WEIGHT                                    =    0 .90203E+01 
VOLUME                                   =    0.  161  08E  +  03 

00 

-0. 

.59688E- 

■02 

3. 

.08 
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TABLE    XXI 

COMPOSITE    DRIVESHAFT    COPES    STANDARD    OUTPUT:  LOGNORMAL 

OPTIMIZATION    RESULTS 

OBJECTIVE    FUNCTION 

GLOBAL    LOCATION  3  FUNCTION    VALUE    0.10281E+02 

DESIGN    VARIABLES 

D.    V.       GLOBAL                     LOWER  UPPER 

ID          NO.       VAR.    NO.                   SOUND                          VALUE  30UND 

1  1                  11  0.1C330E-01         0.19474E-01  0.50000E+00 

2  2                 12  0.10300E-01         0.92916E-01  0.50000E+00 

3  2                  13  0.10000E-01          3.92916E-01  0.50000E+00 

4  3                 14  0.10300E-01         3.10000E-01  0.50000S+00 

5  4                 32  0.0                              0.41392S+02  0.90000E+02 

6  4                 33  0.0                            -3.41392S+02  -0.90000E+02 

7  5                    2  0.50000E+00         0. 2Q465E+0 1  0.5000GE+01 


DESIGN    CONSTRAINTS 


DETERMINISTIC 


GLOBAL  LOWER  MEAN                           UPPER 

ID  VAR.    NO.  BOUND  VALUE                          BOUND 

51  -0.85700E-02  0.37986E-05  0.85700E-02 

3  52  -0.17600E-01  0.3  0.47100E-02 

5  53  -0.18400E-01  0.34484S-02  0.184Q0E-01 

7  54  -0.85700E-02  0.18651E-02  0.85700E-02 

9  55  -0.17603E-01  -0.1S613E-02  0.47100E-02 

11  56  -0.18403E-01  0.46804E-03  0.18400E-01 

13  57  -0.85700E-02  -0.20132S-02  0.85700E-02 

15  58  -0.17600E-01  0.20170S-02  0.47100E-02 

17  59  -0.18400E-01  0.51418E-03  0.18400E-01 

19  60  -0.85700E-02  0.12861E-08  0.85700E-02 

21  61  -0.17600S-01  0.37974E-05  0.47100E-02 

23  62  -0.18400E-01  -0 .  ri  0960E-02  0.18400E-01 

25  63  0.31513E+05  0.54975E+05  0.11000E+16 

26  65  0.30000S+03  0.12564E+04  0.11000E+16 

27  66  -0.11000E+16  0.29994E-01  0.50000E-01 
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TABLE    XXII 

COMPOSITE    DRIVESHAFT    COPES    PROBABILISTIC    OUTPUT :     LOGNORMAL 


PROBABILISTIC 


MEAN 

STANDARD 

COEFFICIENT 

ID 

VALUE 

DEVIATION 

OF    VARIATION 

1 

0-  37986E-05 

0.5215E-06 

0.5215E-01 

3 

0.0 

0.0 

0.0 

5 

0. 34484E-02 

0.5257E-03 

0. 1524E+00 

7 

0. 18651E-02 

0.2842E-03 

0. 1524E+00 

9 

-0.  18613E-02 

0.2838E-03 

0. 15252+00 

11 

0. 46804E-03 

0.2815S-03 

0.6015E+00 

13 

-0.  20132E-02 

0.3069E-03 

0. 1524E+00 

15 

0. 20170E-02 

0.3073E-03 

0. 1524E+00 

17 

0. 51418E-03 

0.3062E-03 

0.5956E+00 

19 

0. 12861E-08 

0.3216E-03 

0. 1000E+02 

21 

0. 37974S-05 

0.3215E-03 

0. 1000E+02 

23 

-0. 40960E-02 

0.6256S-03 

0.1527E+00 

25 

0.54975E  +  05 

0.7813E+04 

0. 1421E+00 

26 

0.  12564E+04 

0.10122+03 

0.8058E-01 

27 

0. 29994E-01 

0.5135E-02 

0. 1712E+00 

PROBABILITY    OF    FAILJ Rj 
AT    LOWER    BOUND 


PROBABILITY    OF    FAILURE 
AT    UPPER    BOUND 


ID 

ALLOWED 

CALCULATED 

ALLOWED 

CALCULATED 

1 

0.  1000E-01 

0.0 

0. 1000E-01 

0.0 

3 

0.  1000E-01 

0.0 

0. 1000E-01 

0.0 

5 

0.  1000E-01 

0.0 

0.  1000E-01 

0.0 

7 

0. 1000E-01 

0.0 

0. 1000E-01 

0.0 

9 

0.  1000E-01 

0.0 

0. 1000E-01 

0.0 

11 

0. 1000E-01 

0.0 

0. 1000E-01 

0.0 

13 

0. 1000E-01 

0.0 

0. 1000E-01 

0.0 

15 

0. 10002-01 

0.0 

0. 1000E-01 

0.3333E-04 

17 

0.  1000E-01 

0.0 

0. 1000E-01 

0.0 

19 

0.1000E-01 

0.0 

0. 1000E-01 

0.5313E-02 

21 

0.  1000E-01 

0.0 

0. 1000E-01 

0.9957E-02 

23 

0. 10002-01 

0.0 

0. 1000E-01 

0.0 

25 

0.  1000E-0  1 

0.1000E-01 

0. 1000E+01 

0.0 

26 

0.1000E-01 

0.0 

0. 1000E+01 

0.0 

27 

0.  1000E  +  01 

0.0 

0. 1000E-01 

0.1544E-02 
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TABLE   XXIII 


COHPOSITE    DHIVESHAFT    DESIGN    SUMMARY:    LOGNORMAL 


COMPOSITE    DRIVESKAFT    OUTPUT 


NUMBER    OF    PLYS  =  % 

NUMBER    OF    MATERIAL    TYPES      =  1 

NUMBER    OF    LOAD    CONDITIONS    =  1 

ECCENTRICITY  =  0.0 


DIMENSIONS 

PLY         THICKNESS  PERCENT         DIAMETER 
INSIDE    0.20465E+01 

1  0.  19474E-01  9.04         0.20855E+01 

2  C. 929162-01  43.16         0.22713E+01 

3  0.92916E-01  43.16         0.24571E+01 

4  0. 10000E-01  4.54          0.24771E+01 


THETA 

0.0 

0.41392E+02 

-0.41392E+02 

0.90000E+02 


STIFFNESS 
AE    =    0.13163E+08 
EI    =    0.82517E+07 
GJ    =    0.95289E+07 


L-C 
1 


LOADS: 
0.31513E+05 


M 
0.0 


0.5000C0E+02    3 


H.P. 
15000E+03 


RPM 
30000E+03 


LOAD 
PLY 

1  0. 

2  0. 

3  -0. 

4  0. 


CONDITION 

E  PL 
37986E-05 
18651E-02 
20132S-02 
12861E-08 


1 

S    F 

100. 00 

4.59 

4.26 

100.00 


EPT 
0.0 

-0.18613E-02 
0.20 170E-02 
0.37974E-05 


S.F. 

130.00 
9.46 
2.34 

100.  00 


EPLT  S.F. 

0.34484E-C2  5*34 

0.46804E-03  39.31 

0.514182-03  35.79 

-0.40960E-02  4.49 


CRITICAL    SPEED  =    0.12564E+04 

MAXIMUM    DEFLECTION  =    0.29994E-01 

WEIGHT  =    0  .  1028 1E+02 

VOLUME  =    0 .  18359E+03 


69 


TABLE   XXIV 
COMPOSITE    DRIVESHAFT    COPES    OUTPUT:    DETERMINISTIC 

OPTIMIZATION    RESULTS 

OBJECTIVE    FUNCTION 

GLOBAL    LOCATION  3  FUNCTION    VALUE    0.  1023 1E  +  02 

DESIGN    VARIABLES 

D.    V.       GLOBAL                    LOWER  U^PER 

ID          NO.       VAR.    NO.                  BOUND  VALUE  EOUND 

1  1                  11                  0.10000E-01  0.10000E-01  0.50000E+00 

2  2                 12                 0.10000E-01  0.10697E+00  0.50000E+00 

3  2                 13                 0.10300E-01  0.10697E+00  0.50000E+00 
a            3                 14                 0.10000E-01  0.10030E-01  0.50000E+00 

5  4                 32                 0.0  0.30524E+02  0.90000E+02 

6  4                 33                 0.0  -0.30524E+02  -0.90000E+02 

7  5                   2                 0.50000E+00  0.18376E+01  0.50000E+01 


DESIGN    CONSTRAINTS 


DETERSI  NISTIC 


GLOBAL  LOWER  MEAN                            UPPER 

ID  VAR.    NO.  BOUND  VALUE                          BOUND 

1  51  -0.42350E-02  0.26961E-05  0.42850E-02 

3  52  -0.88000S-02  0.3  3.23550E-02 

5  53  -0.92030E-02  0.43636E-02  0.92000E-02 

7  54  -0.42850E-02  0.21334E-02  0.42850E-02 

9  55  -0.83000E-02  -0.21307E-02  0.23550E-02 

11  56  -0.92000E-02  0.23560E-02  0.92000E-02 

13  57  -0.42850E-02  -0.23496E-02  0.42850E-02 

15  58  -0.88000E-02  0.23522S-02  0.23550E-02 

17  59  -0.92000E-02  0.26043E-02  0.92000E-02 

19  60  -0.42850E-02  0.17024E-08  0.42350E-02 

21  61  -0.880C0E-02  0.26944E-05  0.23550E-02 

23  62  -0.92000E-02  -0.54219S-02  0.92000E-02 

25  63  0.63025S+05  0.54108E+05  0.11000E+16 

26  65  0.60000E+03  0.13824S+04  0.11000E+16 

27  66  -0.11000E+16  0.24515E-01  0.25000E-01 
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TABLE   XXV 

COMPOSITE    DRIVESHAFT    DESIGN    SUMMARY:    DETERMINISTIC 

COMPOSITE    DRIVESHAFT    OUTPUT 

NUMBER    OF    PLYS  =  % 

NUMBER    OF    MATERIAL    TYPES  =  1 

NUMBER    OF    LOAD    CONDITIONS  =  1 

ECCENTRICITY  =  3.0 

DIMENSIONS 

PLY         THICKNESS         PERCENT         DIAMETER  THETA 
INSIDE    0.18376E+01 

1  0-10000E-01            4.27         0.18576E+01  0.0 

2  0. 10697E+00  45.73         0.20715E+01  0.33524E+02 

3  0. 10697E+00  45.73          0.  22855 E+ 01  -0 . 30 524E+02 

4  0.  10000E-01            4.27         0.23055E+01  0.93000E+02 

STIFFNESS 
AS    =    3.18546E+08 
EI    =    0.99273E+07 
GJ    =    0.66999E+07 


LOADS : 

T     f  T*  M  "P  HO  PPM 

1    *0.31513E+05         0.0      0.500000E+02    3.15000E+03    0.30000E+03 


LOAD    CONDITION  1 

PLY                 EPL  S    F                    ^PT                 S    *•"  PPLT  S    F 

1  0.26961E-05  100.00      0.0                         130 1  00  0.43586E-02  4.*2*1 

2  0.21334E-02  4.02    -0.21307E-02         8.26  0.23560E-02  7.81 

3  -0. 23496E-02  3.65      0.23522E-02         2.00  0.26043E-02  7.07 

4  0.17024E-08  100.00       0.26944E-05    130.00  -0.54219E-02  3.39 

CRITICAL    SPEED  =  0.13824E+04 

MAXIMUM    DEFLECTION  =  3.24515E-01 

WEIGHT  =  0  .  1023 1E*02 

VOLUME  =  0  .  13270S+03 
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VII.   CONCLUSIONS 

Numerical  optimization  using  probabilistic  design  tech- 
niques provides  an  effective  method  for  designing  structures 
based  on  allowed  component  reliability. 

Existing  ANALIZ  subroutines  can  be  easily  modified  to 
perform  probabilistic  optimization  and  as  demonstrated, 
complex  designs  can  be  accomplished. 

The  major  drawback  to  the  program  developed  is  the  large 
increase  in  function  evaluations  as  compared  to  the  standard 
deterministic  optimization.   This  disadvantage  can  be 
partially  offset  by  reducing  the  frequency  in  which  the 
constrained  variable  gradients  are  calculated. 
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VIII.   RECOMMENDATIONS 

A  method  of  scaling  the  variables  should  be  investigated 
to  reduce  ill-conditioning  in  the  optimization  process. 
This  is  a  possible  area  to  pursue  in  order  to  alleviate 
problems  caused  by  very  small  coefficients  of  variation. 

Further  efforts  should  be  undertaken  to  accomplish  the 
reduction  of  function  evaluations  while  maintaining  suitable 
optimization  results.   The  altering  of  the  frequency  in  which 
constrained  variable  gradients  are  calculated  appears  to  be 
the  most  promising  approach. 

In  order  to  reduce  function  evaluations  it  is  recommended 
that  COPES  be  modified  to  allow  the  user  to  supply  mathemati- 
cal expressions  for  the  partial  derivatives  of  the  constrained 
variables  via  the  ANALIZ  subroutine.   These  could  be  used 
instead  of  the  finite  different  gradients  calculated  by 
COPE19/COPE20 .   Modifications  to  allow  the  user  to  supply 
precalculated  gradients  would  also  be  helpful  to  this  end. 
The  sophisticated  user  can  accomplish  this  now  by  using  CONMIN 
directly  in  conjunction  with  C0PE19 ,  C0PE21,  and  COPE22. 

It  is  further  recommended  that  additional  work  be  under- 
taken to  modify  COPES  to  allow  for  the  use  of  the  following 
options  with  probabilistic  design:   A)   two  variable  function 
space;   B)   sensitivity  analysis;   C)   optimum  sensitivity. 
An  additional  worthwhile  modification  would  be  to  provide 
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for  the  calculation  in  C0PE19  and  subsequent  output  of  the 
standard  deviation  of  the  objective  function  mean  value. 

This  investigation  dealt  only  with  component  reliability 
in  the  design  process.   Work  should  be  undertaken  to  ascer- 
tain the  potential  for  the  inclusion  of  system  reliability 
considerations  in  the  optimization  process. 
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APPENDIX  A 
COPES  MANUAL  ADDENDUM 

This  appendix  is  intended  as  an  addendum  to  the  COPES 
Manual,  [Ref .  1] . 

I.   INTRODUCTION 

The  purpose  of  this  document  is  to  provide  user  instruc- 
tions for  performing  numerical  optimization  using  probabilis- 
tic design  techniques.   Four  subroutines  have  been  included 
in  COPES  (Control  Program  for  Engineering  Synthesis)  to 
provide  this  additional  design  tool. 

This  discussion  describes  the  capabilities  of  probabilis- 
tic optimization  using  the  COPES/CONMIN  program.   A  simple 
design  example  is  first  presented  to  demonstrate  the  program 
capabilities.   Guidelines  are  given  for  writing  analysis 
codes.   The  data  organization  is  outlined  and  sample  data  is 
presented. 

This  publication  was  written  to  serve  as  an  addendum  to 
the  COPES  Manual;  "COPES — A  Fortran  Program  for  Engineering 
Synthesis/'  L.  E.  Madsen  and  G.  N.  Vanderplaats ,  NPS69-81- 
00  3,  Naval  Postgraduate  School,  Monterey,  California,  March 
1982. 
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II.   DESIGN  EXAMPLE 

It  is  required  to  design  the  cantilevered  beam  shown  in 
Figure  1.   The  objective  is  to  find  the  minimum  mean  volume 
of  material  which  will  support  the  concentrated  load  and 
maintain  an  allowed  probability  of  failure  of  10%. 


L  =  200  IN— O)  -t>)  B  |<5" 


That  is , 


Figure  1.   CANTILEVERED  BEAM 


Minimize  volume   =   B  *  H  *  L 


(1) 


The  mean  bending  stress  in  the  beam  must  not  exceed  its 

limit  which  is  20,000  psi  which  has  a  coefficient  of  variation 

of  0.0  7; 


BSTRES   = 


Mc 
I 


6  PL 


BH 


2   - 


<   20,000 


(2) 
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The  mean  shear  stress  must  not  exceed  its  limit  of  10,000 
psi  which  has  a  coefficient  of  variation  of  0.05; 


SHRSTR   -   H  -   Sff  -   10'000  (3) 


and  the  deflection  under  the  load  must  not  exceed  one  inch, 
coefficient  of  variation  equals  0.02; 


3         3 
pt,       4pt 
DELTA   =   §|j.  =   i"    <   1.0  (4) 

EBn 


Additionally,  geometric  limits  are  imposed  on  the  mean 
dimensions  so  that; 

0.5   <   B  <_     5.0                       (5) 

1.0   <   H  <   20.0                       (6) 


|   <   10.0  (7) 


The  manufacturing  procedure  is  such  that  all  dimensions  have 
a  coefficient  of  variation  of  0.01.   The  mean  value  of  the 
dead  load  is  10,000  lb.  with  a  coefficient  of  variation  of 
0.07.   The  Young's  modulus  mean  value  is  30.E+0  6  psi  with  a 
0.0  6  coefficient  of  variation. 

The  ANALIZ  subroutine  on  page  11  of  the  COPES  manual 
is  used  for  the  analysis.   Only  one  modification  is  needed. 
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The  variable  P,  for  the  load,  is  added  to  the  global  common 
block. 

The  COPES  data  used  for  the  standard  optimization,  pages 
60  and  61  of  the  COPES  manual  is  modified  to  perform  proba- 
bilistic optimization.  NCALC  in  data  block  B  is  changed  to 
7  to  indicate  probabilistic  optimization  is  to  be  performed 
and  IPROB  equals  1  is  entered  in  column  eight  so  the  normal 
distribution  model  will  be  used.  Three  additional  data 
blocks  are  required,  Ul,  U2  and  U3. 

DATA  BLOCK  Ul : 

There  are  five  probability  variables  which  make  up  the 
constraint  equations  2,  3,  4  and  7.   They  are  B,  H,  E,  AL, 
and  P;  so  IVAR  equals  5.   It  is  desired  that  gradients  of 
the  constrained  variables  only  be  calculated  at  the  beginning 
of  each  CONMIN  iteration,  therefore  LGRAD  equals  1  is  entered 
in  column  two. 

$  DATA  BLOCK  Ul 
5,1 

DATA  BLOCK  U2 : 

The  global  location  of  the  probability  variables  and 
their  corresponding  coefficients  of  variation  are  entered 
in  this  block. 

$  DATA  BLOCK  U2 

1,.01  width,  3 

2,. 01         height,  H 
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10,. 05  load,  P 

8 / • 0  6  Young ' s  mod ,  E 

9,  .01  length,  AL 

DATA  BLOCK  U3: 

This  block  contains  the  allowed  probability  of  failure 
and  coefficient  of  variation  at  both  the  upper  and  lower 
bound  of  each  constrained  variable.   All  four  constrained 
variables  have  upper  bounds,  but  are  unbounded  on  the  lower 
end. 

$  DATA  BLOCK  U3 

1.,.07,.1,.07  BSTRES 

1.,.05,.1,.05  SHRSTR 

1.,.02,.1,.02  DELTA 

1.  , .01, .1, .01  H/B 

The  allowed  probability  of  failure  at  the  upper  bound  and 
the  coefficient  of  variation  are  listed  in  columns  three 
and  four.   Because  there  are  no  lower  bounds  the  resulting 
probability  of  failure  will  be  zero.   These  first  two  columns 
will  be  ignored  during  probability  of  failure  calculations 
in  COPES. 

Coupling  the  ANALIZ  subroutine  to  COPES  yields  the 
following  optimum  design: 
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CANT I LEVERED   BEAM 

AL  =  200.0 

P  =  0.10000E+05 

E  =  0.30000E+08 

B  =  1.92 

H  =  18.75 

VOL  =  7261.54 

BSTRES  =  0.17802E+05 

SHRSTR  =  0.41727E+0  3 

DELTA  =  0.84382E+00 

H/B  =  9.78 

The  COPES  output  for  this  example  is  provided  in  Figure  2. 

III.   PROGRAMMING  GUIDELINES 

The  programming  guidelines  presented  in  Chapter  II  of 
the  COPES  manual  all  apply  in  performing  probabilistic 
optimization.   One  additional  rule  must  be  followed  when 
writing  the  analysis  subroutine:   All  variables  in  the 
constraint  equations  which  have  probabilistic  distributions 
must  be  included  in  the  global  common  block  of  the  analysis 
subroutine. 

IV.   DATA  BLOCKS 

Input  instructions  for  the  probabilistic  optimization 
data  blocks  are  presented  in  the  following  five  pages. 
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V.   NOTES  TO  THE  USER 

1.  The  probabilistic  optimization  feature  is  not 
compatible  with  the  following  COPES  options: 

A.  sensitivity  analysis 

B.  two  variable  function  space 

C.  optimum  sensitivity 

D.  approximate  optimization 

2.  It  is  not  recommended  that  one  attempt  to  design  for 
allowed  probabilities  of  failure  of  less  than  0.00001  or 
greater  than  .99999  due  to  modifications  made  to  "widen"  the 
cumulative  normal  distribution  curve. 

3.  The  optimization  process  is  very  sensitive  to  the 
coefficients  of  variation  of  the  probability  variables  which 
make  up  the  constraint  equations.   Coefficients  of  variation 
of  less  than  0.0001  may  result  in  erroneous  results  if  the 
coefficients  of  variation  of  some  of  the  constraint  limits 
are  also  less  than  0.0001. 

4.  In  some  cases  the  initial  design  will  be  such  that 
the  optimization  process  terminates  unsatisfactorily  due  to 
the  calculation  of  zero  constraint  gradients.   Termination 
will  occur  either  due  to  an  inability  to  achieve  a  feasible 
design  or  CONMIN  will  reach  the  maximum  number  of  iterations 
no  matter  how  large  ITMAX  is  made.   This  probably  will  be 
due  to  the  fact  that  the  probabilities  for  failure  calcu- 
lated are  either  1.0  or  0.0  and  the  finite  difference  steps 
are  insufficient  to  reach  the  region  along  the  cumulative 
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normal  distribution  curve  between  Pf  =  1.0  and  Pf  =  0.0. 
This  could  be  caused  by  an  initial  design  which  is  a  great 
distance  from  the  optimum  or  by  very  small  coefficients  of 
variation.   If  this  happens,  perform  a  deterministic  optimi- 
zation to  achieve  a  design  in  the  vicinity  in  the  design 
space  desired  and  use  this  result  for  the  initial  design  in 
probabilistic  optimization. 

5.   In  order  to  determine  probabilities  of  failure,  COPES 
calculates  the  gradients  of  the  constrained  variables  with 
respect  to  the  probability  variables  which  make  up  the 
constraint  equations.   This  is  done  by  finite  difference 
steps  resulting  in  numerous  calls  to  the  ANALIZ  subroutine. 
In  many  problems  these  gradients  remain  essentially  the  same 
so  the  user  is  given  the  ability  to  determine  the  frequency 
in  which  these  gradients  are  calculated.   This  provides  for 
a  significant  reduction  in  the  number  of  function  evaluations 
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